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WAVELET FILTERING 
Technical Field of the Invention 

The present invention relates generally to wavelet filters and, in particular, to 
wavelet filtering of digital images, 
5 Background 

The field of digital data compression and in particular digital image compression 
has atti'acted great interest for some time. 

In the field of digital image compression, many different techniques have been 
utilised, Li particular, one popular technique is the JPEG standard, which utilises the 
10 discrete cosine transform to tratisform standard size blocks of an image into 
corresponding cosine components. In this respect, the higher frequency cosine 
components are heavily quantised so as to assist in obtaining substantial compression. 
The heavy quantisation is an example of a lossy technique of image compression. The 
JPEG standard also provides for the subsequent lossless compression of the tratisform 
15 coefficients. 

Recentlyj the field of wavelet transforms has gained great attention as an alternative 
form of data compression. The wavelet transform has been found to be highly suitable in 
representing data having discontinuities such as sharp edges. Such discontinuities are 
often present in image data or the like, 

20 Typically, data compression using wavelet techniques is a Uvo step process. It 

comprises, firstly, a transform phase, during which the wavelet transform of the data set is 
calculated, and secondly a subsequent coding stage during which the resultant data set 
from the transform operation is separated into segments which are then coded using a 
specific coder. In decompression, the reverse occurs, with coded blocks being first 

25 decoded, and subsequently the inverse wavelet transform being applied to generate the 
final decompressed output. 

Typically, the wavelet transform operation utihses a number of wavelet filters. 
These filters are typically digital filters, such as finite impulse response (FIR filters) 
which ideally are linear and require a small number of taps. These filters process the 

30 digital images in bands and typically require a large external store for storing overlapping 
pixels, coefficients, and/or intermediate results for processing the next band. The 
bandwidth to the external store is typically limited, which restricts the rate of reading mid 
writing of overlapping pixels, coefficientSj and/or intermediate results to the store. 
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Disclosure of the Invention 

It is an object of the present invention to substantially overcome, or at least 
ameliorate, one or more disadvantages of existing wavelet filters. 

According to a first aspect of the invention, there is provided a method of wavelet 
5 filtering a digital image, the digital image comprising a plurality of pixels arranged in one 
or more bands of pixels, the method performing the following steps for each one of said 
one or more bands of said pixels: retrieving a plurality of previous partial results from a 
first or second storage; inputting a current group of adjacent said pixels; computing 
coefficients utilising said current group of adjacent pixels and said plurality of previous 

10 partial results; outputting said coefficients; and storing a plurality of current partial results 
in said first storage if said current group is not the last group in the band, in said second 
storage if said group is the last group in the band. 

According to a second aspect of the invention^ there is provided an apparatus for 
wavelet filtering a digital image, the digital image comprising a plurality of pixels 

15 arranged in one or more bands of pixels, the apparatus comprising: means for retrieving, 
for each one of said one or more bands of pixels, a plurality of previous partial results 
from a first or second storage; means for inputting, for each one of said one or more 
bands of pixels, a current group of adjacent said pixels; means for computing coefficients, 
for each one of said one or more bands of pixels, utilising said current group of adjacent 

20 pixels and said plurality of previous partial results; means for outputting, for each one of 
said one or more bands of pixels, said coefficients; and means for storing a plurality of 
current paitial results in said first storage if said current group is not the last group in the 
band, in said second storage if said group is the last group in the band. 

According to a third aspect of the invention, there is provided a computer readable 

25 medium comprising a computer program for wavelet filtering a digital image, the digital 
image comprising a plurality of pixels arranged in one or more bands of pixels, the 
computer program comprising: code for retrieving, for each one of said one or more 
bands of pixels, a plurality of previous partial results from a first or second storage; code 
for inputting, for each one of said one or more bands of pixels, a current group of adjacent 

30 said pixels; code for computing coefficients, for each one of said one or more bands of 
pixels, utilising said current group of adjacent pixels and said plurality of previous partial 
results; code for outputting, for each one of said one or more bands of pixels, said 
coefficients; and code for storing a plurality of current partial results in said first storage 
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if said current group is not the last group in the band, in said second storage if said group 
is the last group in the band. 

According to a fourth aspect of the invention, there is provided a wavelet filter for 
wavelet filtering a digital image, the digital image comprising a plurahty of pixels 
5 arranged in one or more bands of pixels, the filter comprising: a first storage for storing a 
plurality of previous partial results; a second storage for storing a plurality of previous 
partial results^ a controller for selecting the plurality of previous partial results from said 
first or second storage; a pixel input mechanism for inputting a current group of adjacent 
said pixels; a lifting lattice of multiplier and adder units for computing coefficients 
10 utilising said current group of adjacent pixels and said selected plurality of previous 
partial results; output means for outputting said coefficients; and a controller for storing a 
plurality of current partial results in said first storage if said current group is not the last 
group in the band, in said second storage if said group is the last group in the band. 

Brief Description of the Drawings 
15 A number of preferred embodiments of the present invention will now be described 

with reference to the drawings^ in which: 

Fig. 1 is a schematic representation of a K-line 5-3 wavelet filter (with K = 8) for 
filtering a digital image in accordance with a first embodiment; 

Fig. 2 is a schematic representation of the cascading of a wavelet filter of Fig. 1 in 
20 time; and 

Fig, 3 is a flow diagram of a method for use in wavelet filtering in accordance with 
a second embodiment; and 

Fig. 4 is a schematic representation of a general-purpose computer for use in 
implementing the method of Fig. 3. 
25 Detailed Descriptiod including Best Mode 

Where reference is made in any one or more of the accompanying drawings to steps 
and/or features, which have the same reference numerals, those steps and/or features have 
for the purposes of this description the same function(s) or operation(s), unless the 
contrary intention appears. 
30 The principles of the preferred embodiments described herein have general 

applicability to the implementation of K-line wavelet filters. However^ for ease of 
explanation, the embodiments has been described with reference to a K-line 5-3 wavelet 
filter for use in a forward discrete wavelet transform of digital images. However, it will 
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be readily evident that the invention is not limited thereto. For instance, the wavelet filter 
can also be used in K-line 9-7 wavelet filters. Furthennore, it would be apparent to a 
person skilled in the art, that the wavelet filter can be modified for use in an inverse 
discrete wavelet transform. Still further, it would be apparent to a person skilled in the art 
5 that two such wavelet filters can be combined in known manner to provide a 2- 
dimensional discrete wavelet transfomi or modified for use in a 2-dimensional inverse 
discrete wavelet transform. Moreover, the wavelet filter can be used in many different 
applications. For examples of the many different applications of wavelet analysis to 
signals, reference is made to a survey article entitled "Wavelet Analysis" by Bruce et al. 
10 appearing in IEEE Spectrum, October 1996, pages 25 to 26. For a discussion of the 
different applications of wavelets in computer graphics, reference is made to 'Wavelets 
for Computer Graphics", page 5, L Stollinitz et, al. published 1996 by Morgan Kaufiiami 
Publishers, Inc. 

As mentioned above, the embodiments of the invention can be used in a discrate 
15 wavelet filter transform as well as an inverse discrete wavelet transfoim. For sake of 
simplicity^ the tenn wavelet transform used herein or variations thereof, is taken to 
include both a forward and/or an inverse wavelet transform, unless the contrary intention 
appears. Similarly, the term pixel as used herein or variations thereof, is taken to taken to 
refer to the original pixels or wavelet coefficients of a digital image, unless the contrary 
20 intention appears. Similarly, the term digital image or variations thereof is taken mclude 
an original image or a sub-band of that original image, or one or more associated sub- 
bands of wavelet transformed coefficients of that original image. 

Turning now to Fig. 1, there is shown a schematic representation of a K-line 5-3 
wavelet filter (with K - 8) for filtering a digital image in one dimension in accordance 
25 with a first embodiment. The digital image comprises a plurality of pixels arranged in a 
plurality of columns in the vertical direction and a plurality of rows in the horizontal 
direction. > 

The filter 100 processes one horizontal band of the digital image after another. The 
digital image comprises one or more such horizontal bands, with the horizontal bands 
30 comprising a plurality of said rows of said pixels. The filter 100 initially commences 
processing the uppermost band of the image, and then processes the next adjacent band of 
the image, and so on until the lowermost band of the image is processed. Within each 
band, the filter 200 processes one column of pixels after another. The filter 100 first 
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processes the left column of the band, and then the next adjacent column of the band, and 
so on until the right column of the band is reached. For ease of understanding, the pixels 
in each column of a baud from the topmost pixel to the lowermost pixel are represented as 
^0,^1.^2. - -- -^a:-/' 

5 Within each column of a band, the filter 100 processes a first group of adjacent 

pixels (commencing with the topmost group) and then processes a second group of 
adjacent pixels adjacent the first group. 

In this way, the filter 1 00 performs a one dimensional discrete wavelet transform in 
the vertical direction. In another embodiment, the filter 100 can also be adapted to 

10 perform a one dimensional discrete wavelet transform in the horizontal direction by 
processing vertical bands one after another in similar fashion. In still another 
embodiment, the filler can also be adapted to provide a 2-dimensional wavelet transfonn. 
In the latter case, one filler 100 processes one column at a time first, and a second filter 
100 after that collects the output of the first filter 100 and processes one row at a time. 

15 The pixels ( e.g. xo, xj^ X2, xs^ Xs^ Xc^ xj^ xs, xp^ xjo, xn, x/s, x/if of a column 

of a band are fed to the filter 1 00. These pixels are fed in parallel to the wavelet filter 100 
via a number 2K (with K = S) of input lines 101, The 16 input lines 101 are connected to 
eight 2-input multiplexers 103-1, 103-2, 103-3, 103-4, 103-5, 103-6, 103-7, and 103-8, 
The multiplexers can select either the bottom eight lines or top eight lines for input of the 

20 groups of pixels (pco^ x/^ X2, xj, X4, xs, Xo, and ^7) or (xg^ X9, xjq^ xjj^ xn, xj4 and xjs) 
respectively. The multiplexers 103-1, 103-2, 103-3, 103-4, 103-5, 103^6, 103-7, and 103- 
8 are connected to a lifting lattice of multiplier and adder units 102, 104, 106, 108, 110, 
112, 114, and 116, which compute the high pass coefficients and low pass coefficients on 
the top 8 lines (ie. coefficients (xs^ xg^ xjo^ xjj^ xj2, xjj^ xj4 and xji}) or the bottom S lines ( 

25 ie. coefficients {xo, x/, X2, X3, X4^ xs, xc, and x?))* The lifting lattice of multiplier and adder 
units 102, 104, 106, 108, 110, 112, 114, and 116, perfomi the following calculations: 

<^n =^2« +^^(^2.-1 « = 0,1,„.,A:"1 (1) 

and, 

30 = x^,,, + /3{d„ + d,,,, ) n = -1A^.K - 2 (2) 

where d,, are all the high pass coefficients, s„ all the low pass coefficients, a = -0.5, and (5 
= 0.25, 
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For example, the multiplier and adder units 102, 104 106, 108, 110, 112, 114, and 
116 perform the following calculations on the selected input pixels (xo, xj, X2, x^. x^^ xs x^^ 
and xt) to compute the high pass coefficients dn and the low pass coefficients Sni 



5 


unit 102 


^0 = ^0 + ^(^-1 + ^3 ) 


(3) 




unit 104 




(4) 




unit 1 06 


^2 = ^4 + CC{X2 + ^5 ) 


(5) 




unit lOS 




(6) 




unit 110 




(7) 


10 


unit 112 




(8) 




unit 1 14 




(9) 




linit U6 


^2 ^^5^/^{d2 ^d,) 


(10) 



The high pass and low pass coefficients dn and s„ are then passed to other discrete 
J 5 wavelet transform circnits or entropy coder (not shown in the figure). 

The multiplier and adder units 110 and 1 02 also accept inputs from either the partial 
result registers 109 or external buffer 107, depending on the control on the 
multiplexer 105. Namely, the multipKer and adder units 110 and 102 accept input from 
the external buffer ] 07 when computing the first 8 lines and accepts input from the partial 
20 result registers 109 for the last 8 lines. 

A controller 1 1 1 controls the storing of a last input pixel and a last high pass 
coefficient (e.g. xy and ) to the partial result register 109 or external buffer 107 
depending on which lines are being processed. Namely, the last input pixel and the last 
high pass coefficient of the first 8 lines are stored in the partial result register 109. 
25 However, when the lifting lattice of multiplier and adder units is computing coefficients 
for the last S lines, the last input pixel and last high pass coefficient are stored in the 
external buffer 107. The last input pixel and last high pass coefficient are required for 
computation of the low and high pass coefficients of the next adjacent group of pixels 
within the column, 

30 For example, in the case where the filter 100 is currently processing the pixels {xs. 

X9. xio, x,i. x,2. xn. xi4 , and xis), the multiplier and adder ujiits 1 10 and 102 accept as input 
the last pixel and coefficient x? and di from the partial result registers 109. which were 
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previously stored in the partial result register 109 during the previous computation of the 
pixels (xo. xj, X2, xs^ X4, xs, xe, and xj) of the same column and band. In the case where the 
filter 100 is currently processing the pixels (xo to X7% the multiplier and adder units 110 
and 102 accept as input the last pixel and coefficient Xi$ and from the external 
5 buffer 107. These inputs xis and d? were previously stored in the external buffer 107 
during the previous computation of the pixels {xsto xjs) of the same coluni^n but in the last 
previous band. In the case where there are no previous pixels (e.g. topmost band) the 
boundary is synunelrically reflected- 

Figure 2 illustrates the cascading of the wavelet filter of Fig. 1 in time. In cycle 0, 

10 for the jfirst set of inputs the wavelet filter 100 either symmetrically reflects the pixels 
near the boundary or reads tlie required partial results firom external buffer 107. The 
filter 100 then stores the partial results into the partial results register 109, and outputs the 
low and high pass coefficients. In cycle 1, the lifting lattice is 'moved' to the next set of 
input lineSj reads input from partial results register 109, and outputs the next low and high 

15 pass coefficients. The partial results are stored in the partial results register 109 and 
external buffer 107. 

The first embodiment has been described with respect to a filter 100 having two 
groups of eight input lines 101 and eight two-input multiplexers 103-1 to 103-S, It would 
be apparent to a person skilled in the art^ that the mnnber of groups can be increased (e.g. 

20 N = 4 groups) with a corresponding modification of the eight multiplexers, (ie. eight N- 
input multiplexers). The external buffer 107 is a slow memory storage and requires at 
least a memory size equivalent to the size of an input pixel and a high pass coefficient 
times the nuniber of columns. The filter 1 00 stores in external memory the last input 
pixel and the last high pass coefficient xjs and dr of each column of each band for fi.irther 

25 processing in the next band. On the other hand, the partial results register is a fast local 
storage and requires only a memory size of an input pixel and high pass coefficient. 

The effective value of K' (the width of the band) can be increased by increasing the 
number of groups of pixels (eg. two groups, K* = 2xK ^ 16) or (four groups K' = 4xK 
^32)). In the present embodiments, this has the advantage of reducing the number of 

30 accesses to the slow external memory. Namely, where there are multiple groups, the last 
pixel and last high pass coefficient can be stored in the fast local memory, namely the 
partial results register 109, rather than the slow external memory 107, thus enabling the 
wavelet transform to be performed more efficiently. Of course, the last pixel and last 
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high pass coefficient of the last group is required to be stored in external memory. 
Furthermore, the amount of hardware required to implement the filter is not significantly 
increased. Still further, the filter has the capability of filtering with different values of K\ 
In an alternate embodiment, the multiplexers 103-1, 103-2, 103-3, 103-4, 103-5, 
5 103-6, 103-7^ and 103-8 may be dispensed with. In this case, the multiplexing is done by 
pixel fetching mechanism implicitly. For instance, the pixels are fed to the lifting lattice 
from a frame buffer by an address generator. 

Turning now to Fig. 3, there is shown a method for use in wavelet filtering in one 
dimension in accordance with a second embodiment. The method comprises in part a 

10 sub-procedure 300 for use in the wavelet filtering of a digital image. As mentioned 
previously, the digital image comprises a plurality of pixels arranged in a plurality of 
columns in the vertical direction and a plurahty of rows in the horizontal direction. The 
digital image comprises one or more horizontal bands, each of which comprise a plurality 
of said rows of pixels. The method processes one horizontal band of the digital image 

15 after another. The method initially commences processing the uppermost band of the 
image, and then processes the next adjacent band of the image, and so on until the 
lowermost band of the image is processed. Within each band, the method processes one 
column of pixels after another. The method first processes the left column of the band, 
and then the next adjacent column of the band, and so on until the right column of the 

20 baad is reached For ease of understanding, the pixels in each column of a band from the 
topmost pixel to the lowermost pixel are represented as xo,xj^ Xj, . . , xk'-i- 

The method of wavelet filtering calls the sub-procedure 300 for processing each 
column of a band. Within each column of a band, the sub-procedure 300 processes in 
tura a first group of K adjacent pixels^ a second group of K adjacent pixels adjacent the 

25 first group , a third group of 5C adjacent pixels adjacent the second group, and so on to the 
last group of pixels. The top and bottom pixel of the first and last group of pixels 
respectively, being located adjacent the boundary of the band. 

In this way, the method of wavelet filtering performs a one dimensional discrete 
wavelet transform in the vertical direction. In another embodiment, the method can also 

30 be adapted to perform a one dimensional discrete wavelet transfomi in the horizontal 
direction by processing vertical bands one after another in similar fashion. In still another 
embodiment, the method can also be adapted to provide a 2-dimensional wavelet 
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transform. In the latter case, the method processes one colunm at a time first, and the 
method after that collects the output and processes one row at a time, 

The sub-procedure commences at step 301 and proceeds to step 302 where a 
counter M is set to one. As mentioned previously, the sub-procedure processes the pixels 
within a column of a band. The sub-procedure continues to a loop 306, 308, 310, 312, 
314, and 316, where groups of a number K adjacent pixels are retrieved and processed 
one after another. A first group of K adjacent pixels (nearest the top boundary) is 
processed during the first pass of the loop 304 to 313, a next group of K adjacent pixels 
(adjacent the first group) is processed during the second pass of the Joop 304 to 313, and 
so on The processing of one group of pixels after another is achieved by incrementing 
the counter M. The sub-procedure terminates once K' pixels have been processed, when 
M.K - K\ Preferably, K = 4 and 4 or S. 

After the counter M has been set to one, the sub-procedure proceeds to step 304. hi 
step 304, the method reads the corresponding previous partial results from an external 
buffer. These previous partial results having being stored in the external buffer during the 
processing of the sub-procedure of the same column in the last previous band of pixels. 
In the case where there are no previous pixels (e.g. topmost band) the boundary is 
symmetrically reflected. 

After step 304, the method proceeds to loop 304 to 316. This loop will be 
described with respect to the general case where M = M'. During step 306, the sub- 
procedure reads the M'th group of K adjacent pixels and then computes the low and high 
pass coefficients in step 308 in accordance with the equations (1) and (2) described 
above. In computing these high and low pass coefficients the method takes as input the 
M'th group of K adjacent pixels. It also accepts input from either the external buffer or a 
local register which stores the last pixel and the last high pass coefficient. The step 308 
outputs to other discrete wavelet transform or entropy encoding procedures (not shown in 
the figure) . 

For example, in the case where the sub-procedure is currently processing the pixels 
(X8, xg^ xjo. x/j^ xj2, x;j^ Xj4 , and xjs), the sub-procedure step 308 accepts as input the last 
pixe] and high pass coefficient x? and ds from the partial result registers, which were 
previously stored in the partial result register during the previous computation of the 
pixels {xo^ XI ^ X2, X4, x^^ xts^ and x?) of the same cohunn and band, hi the case where the 
sub-procedure is currently processing the pixels (xo to X7X the sub-procedure step 308 
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accepts as input the last pixel and high pass coefficient from the external buffer. These 
latter inputs were previously stored in the external buffer during the previous computation 
of the last group of pixels of the same column but in the last previous band. In the case 
where there are no previous pixels (e.g, topmost baud) the boundary is symmetrically 
5 reflected- 

Ater step 308, the method proceeds to a decision block 310, where a check is made 
whether M.K = K\ If the decision block returns false, then the method proceeds to 
step 312j where the last pixel and last high pass coefficient are stored in the local partial 
results register for use in the computation 308 of the next group of pixels. 

10 In the next step 314, the counter M is incremented and the last pixel and the last 

high pass coefficient are read 316, The method then returns to step 306 for the next pass 
of the loop 306 to 316, In the event, the decision block 310 returns true, then the last 
pixel and the last high pass coefficient are stored in the external buffer and the sub- 
procedure finishes for the current column of the band, 

15 As can be seen, the loop can be repeated as many times as desired to increase the 

value of K\ 

The aforementioned preferred method comprise a particular control flow. There are 
many other variants of the preferred method which use different control flows without 
departing the spirit or scope of the invention. Furthermore one or more of the steps of the 

20 preferred method may be performed in parallel rather sequential. 

The method of Fig. 3 is preferably practiced in dedicated hardware, such as shown 
in Fig, 1. However, the method of Fig, 3 can also be practiced using a conventional 
general-purpose computer system 400, such as that shown in Fig. 4 wherein the processes 
of Fig. 3 may be implemented as software code, such as an application program executing 

25 within the computer system 400, In particular, the steps of method of Fig. 3 are effected 
by code instructions in the software that are carried out by the computer. This method 
has particular application in computer systems restricted to large slow external buffers 
and small fast local memories, 

The software may be stored in a computer readable medium, including the storage 

30 devices described below, for example. The software is loaded into the computer from the 
computer readable medium, and then executed by the computer 
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The use of the computer program readable medium together with the software in 
the computer system preferably effects an advantageous apparatus in accordance with the 
embodiments of the invention. 

The computer system 400 comprises a computer module 401, input devices such as 
5 a keyboard 402 and mouse 403, output devices including a printer 41 5 and a display 
device 414. A Modulator-Demodulator (Modem) transceiver device 416 is used by the 
computer module 401 for communicating to and from a communications network 420, for 
example connectable via a telephone line 421 or other functional medium. The 
modem 416 can be used to obtain access to the Internet, and other network systems, such 

10 as a Local Area Network (LAN) or a Wide Area Network (WAN). 

The computer module 401 typically includes at least one processor unit 405, a 
memory unit 406, for example formed from semiconductor random access memory 
(RAM) and read only memory (ROM), input/output (I/O) interfaces including a video 
interface 407^ and an I/O interface 413 for the keyboard 402 and mouse 403 and 

15 optionally a joystick (not illustrated), and an interface 408 for the modem 416. A storage 
device 409 is provided and typically includes a hard disk drive 410 and a floppy disk 
drive 411. A magnetic tape drive (not illustrated) may also be used. A CD-ROM 
drive 412 is typically provided as a non- volatile source of data. The components 405 
to 413 of the computer module 401^ typically communicate via an interconnected bus 404 

20 and in a manner which results in a conventional mode of operation of the computer 
system 400 known to those in the relevant art. Examples of computers on which the 
embodiments can be practised include IBM-PC's and compatibles. Sun Sparcstations or 
alike computer systems evolved therefrom. 

Typically, the application program of the preferred embodiment is resident on the 

25 hard disk drive 410 and read and controlled in its execution by the processor 405. 
Intermediate storage of the program and any data fetched from the network 420 may be 
accomplished using the semiconductor memory 406, possibly in concert with the hard 
disk drive 410. In some instances, the application program may be supplied to the user 
encoded on a CD-ROM or floppy disk and read via the corresponding drive 412 or 411, 

30 or alternatively may be read by the user from the network 420 via the modem device 416. 
Still further, the software can also be loaded into the computer system 400 from other 
computer readable medium including magnetic tape> a ROM or integrated circuit, a 
magneto-optical disk, a radio or infra-red transmission channel between the computer 
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inodule 401 and another device, a computer readable card such as a PCMCIA card^ and 
the Internet and Intranets including email transmissions and information recorded on 
websites and the like. The foregoing is merely exemplary of relevant computer readable 
mediums. Other computer readable mediums may be practiced without departing from 
the scope and spirit of the invention, 

The aforementioned embodiment shown in Fig. 1 can be modified for use in a 5/3 
inverse transform wavelet filter. In this modified embodiment, the 5/3 inverse transfonn 
wavelet filter calculates the original pixels from the low pass coefficients and high pass 
coefficients in accordance with the following formulae: 

Pid„.,'^dJ n = OJ,...,K-l (11) 

and, 

^2..i =dn- «(^2« + ^1.^2 ) n = -1,0,.,^: - 2 (12) 
where a = -0,5 and (3 = 0,25, The multiplier and adder units of Fig, 1 will be modified 
accordingly to compute the pixels in accordance with equations (1 1) and (12). This filter 
takes as input the coefficients sq to s^-i and do to d^^-j ^ and also requires as input the 
coefficients d.f and x.2 from the previous band^ which are stored in the partial register 109 
or external buffer 107. 

The aforementioned embodiment of Fig, 1 can also be modified without departing 
fi:om the spirit of the invention to other wavelet filters such as a 9/7 wavelet filters. 

In this modified embodiment, the 9/7 forward transform wavelet filter computes the 
low pass coefficients and high pass coefficients in accordance with the following 
formulae: 

^'n=^zn.i^ a(x^^'hX^^^J n = 0,1,...,K-1 (14) 
^n^X2o^ P(d'n^,'^d'J n = 0,1,„„K-1 (15) 

<in =dn-^r( S„ •^Sn.i ) n= -^0,,..,^- 2 (16) 

3, =3\+S(d,_, +dj n = '1A.:.K - 2 (17) 
where d„ are all the high pass coefficients, s„ all the low pass coefficients, d'n and s « are 
intermediate values, a = -1.5861, P = -0.052980, y = 0.88291, and 5= 0.44351. The 
multiplier and adder units of Fig. 1 will be modified accordingly to compute the 
coefficients in accordance with equations (14) to (17). This filter takes as input pixels xj 
to x^ and also requires as input the pixel xo, the intermediate values d'.j and s\j, and the 
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high pass coefficient d,z from the previous band, which are stored in the partial register 
109 or external buffer 107, 

In another modified embodiment, the 9/7 inverse transform wavelet filter computes 
the original pixels from the low pass coefficients and high pass coefficients in accordance 
with the following formulae: 

s;=s,- ^rc/..,+tfJ n^0.t,.,,K-1 (18) 
=d, - r(s, ) n ^-1A-.K - 2 (19) 
x,,=s\,- p(d,,,^d\) n=-1A^,^.K-2 (20) 

Xz..i =dn -<^( ^^20^2) " = -2,-7,.,.,/< -3 (21) 
where dn are all the high pass coefficients, Sn all the low pass coefficients, d\ mds\ are 
intermediate values, a = -L5861, (5 - -0.052980, y = 0.88291, and 5= 0.44351. The 
multiplier and adder units of Fig. 1 will be modified accordingly to compute the pixels in 
accordance with equations (1 8) to (21), This filter takes as input coefficients sq to Sic^j and 
do to dfc-i and also requires as input the pixel x^, the intermediate values d\2 and 5 and 
the high pass coefficient d.j from the previous band which are stored in the partial 
register 109 or external buffer 107. 

The aforementioned embodiment shown in Fig. 3 can also be modified in similar 
fashion to other wavelet filters, such as a 5/3 inverse transform wavelet filter or 9/7 
wavelet filters. 

Industrial Applicability 

It is apparent from the above that the embodiment(s) of the invention are applicable 
to the computer and data processing industries. In particular, the embodiments have 
application in digital data compression of digital images and videos. 

The foregoing describes only one embodiment/some embodiments of the present 
invention, and modifications and/or changes can be made thereto without departing from 
the scope and spirit of the invention, the embodiment(s) being illustrative and not 
restrictive. 
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